Hall plateau diagram for the Hofstadter butterfly energy spectrum 
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We extensively study the localization and the quantum Hall effect in the Hofstadter butterfly, 
which emerges in a two-dimensional electron system with a weak two-dimensional periodic potential. 
We numerically calculate the Hall conductivity and the localization length for finite systems with 
the disorder in general magnetic fields, and estimate the energies of the extended levels in an 
infinite system. We obtain the Hall plateau diagram on the whole region of the Hofstadter butterfiy, 
and propose a theory for the evolution of the plateau structure with increasing disorder. There 
we show that a subband with the Hall conductivity ne^ /h has \n\ separated bunches of extended 
levels, at least for an integer n < 2. We also find that the clusters of the subbands with identical 
Hall conductivity, which repeatedly appear in the Hofstadter butterfly, have a similar localization 
property. 



I. INTRODUCTION 

A two dimensional (2D) electron system with a 2D pe- 
riodic potential is expected to exhibit an intricate energy 
spectrum in a strong magnetic field, which is called the 
Hofstadter butterfly. When the magnetic length is of the 
order of the lattice constant, the interplay of the Landau 
quantization and Bragg reflection yields a fractal-like se- 
ries of the energy gaps, which depends sensitivelyon the 
number of magnetic flux quanta per a unit cell, 0] 



^ h/e' 

where B is the amplitude of the constant magnetic field 
and a is the lattice constant. Moreover it is shown 
that each subband carries a quantized Hall conductivity, 
which varies with the energy gaps in a nontrivial manner. 
The recent developments in experimental techniques 
make it possible to fabricate a two-dimensional super- 
lattice on a 2D electron gas. Evidence of the fractal 
energy spectrum was found in a superlattice patterned 
on GaAs/AlGaAs heterostructures, where the Hall 
conductivity changes nonmonotonically as the Fermi en- 
ergy transfers from one gap to another. 

From the theoretical side, it is intriguing to consider 
the magnetotransport in this intricate energy spectrum. 
Broadening of the density of states and the con- 
ductivity |5i were investigated for the Hofstadter butter- 
fly within the self-consistent Born approximation. For 
the localization regime, we expect that, by analogy with 
unmodulated 2D systems, the extended levels appear 
only at certain energies in the Hofstadter spectrum, and 
Cxy (E) will turn into a series of Hall plateaus separated 
by those energies. 

The evolution of the extended states as a function of 
the disorder was proposed for several flux states in the 
Hofstadter butterfly. |3j S Pi A finite-size scaling analy- 
sis was performed for a 2D system modulated by a weak 
periodic potential and it was found that the modulation 
does not change the critical exponent at the center of 



the Landau level. [13, 113 Previously we have numeri- 
cally calculated the Hall conductivity in weakly modu- 
lated 2D systems at several fiuxes and studied the effect 
of the localization on the quantum Hall effect in the Hof- 
stadter butterfly. (12] There we determined the energies 
of the extended levels from the scaling behavior of the 
Hall conductivity in finite systems. 

In this paper, following the line of the previous work, 
we present an extensive study of the quantum Hall ef- 
fect in a single Landau level in the presence of the two- 
dimensional periodic potential. We numerically calculate 
the Hall conductivity and the localization length in gen- 
eral magnetic fluxes, using the exact diagonalization of 
the Hamiltonian for finite systems. We determine the en- 
ergies of the extended levels in an infinite system by an- 
alyzing the size dependence in the Hall conductivity, and 
find that those energies always coincide with the points 
where axy/{—e'^/h) becomes n -I- 1/2 (n is an integer) in 
finite systems. In particular, we resolved two separated 
extended levels in a subband carrying the Hall conduc- 
tivity by 2e^/h. We propose a possible scenario for the 
evolution of the extended levels with increasing disorder 
for several fluxes, which suggests that the Hall plateau 
structure changes through the pair annihilation and/or 
pair creation of the extended levels. We survey the Hall 
conductivity over the whole region of the Hofstadter but- 
terfly, to find that the electronic structure with similar 
Hall conductivity and localization length appears in a 
self-similar manner at different fluxes, reflecting the frac- 
tal property of the ideal spectrum. 



II. FORMULATION 

We first prepare the formulation to describe a Bloch 
electron in magnetic fields. Let us consider a two- 
dimensional system in a uniform magnetic field with a 
periodic potential Vp and a disorder potential Vd, 
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(p + eAf + Vp + Vd^ 



(1) 
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We assume that Vp has a square form 

\/„ = \/ cos — x+V COS — y, 
a a 

and the disorder potential is composed of randomly dis- 
tributed S potentials ±uo with the number per unit area 
Hi , where the amounts of the positive and negative scat- 
terers are taken to be equal. We consider only the lowest 
Landau level (A^ = 0), assuming that the magnetic field 
is strong enough and the mixing of the Landau levels is 
neglected. In the Landau gauge, the basis can be taken 
as 



|0,fc,) 



1 



exp{ikyy) exp 



2/2 



(2) 



with the magnetic length I — ^ h/eB. The matrix ele- 
ments of Vp are then written 

{G,k'y\Vp%ky) = 5k-^,k2Ve-^^cos^ 

+ (4;,fc„+27r/a + 5k'^^ky-2T:/a)Ve~^ . (3) 

In an ideal system {Vd — 0), the wave function can 
then be expanded as 



"^kyir) = ^Cm(fcj,)|0, ky 



271- 



(4) 



with Bloch wave number ky ranging from — vr/a to 7r/a. 
The Schrodinger equation is then reduced to Harper's 
equation, 

Fe--/2'^(c„,+i-|-c„._i) 

+2Ve--'^^ cos ( = Ee^. (5) 



In a rational flux </> = p/q {p,q are coprime integers), 
the equation becomes periodic in m with a period p, 
which corresponds to a distance qa in the center coor- 
dinate along the x axis, so that we have the Bloch con- 
dition Cm+p = e*'^"='''Cm with another Bloch wave num- 
ber kx from — 7r/(ga) to T:/{qa). As a result, we have 
p independent states for each of A; = (kx,ky) in the 
Landau level, so that we can label the wave function as 
V^nfe on a g-folded Brillouin zone with the subband index 
n = 1,2, We can decompose the wave function as 

'0nfc('') = e^'^^^'^^^^y^ Unk{f), where u satisfies the mag- 
netic Bloch condition 



L{x,y + h) = u{x,y) 



u{x + qa,y) 



u{x,y). 



(6) 



Figuren]shows the energy spectrum in the lowest Lan- 
dau level at Vd = plotted against cj). The intricate band 
structure is due to the number of subbands p, which is 
not a continuous function of 0. We can show that p sub- 
bands never overlap so that we always have p — I energy 



gaps inside the Landau level. The total width of the spec- 
trum scales with the factor e^'^/^^ and shrinks 
as the flux becomes smaller. The lower panel shows the 
zoom out of the spectrum covering from = to 10. 
We can see that, on going to a higher field, a series of 
subbands splits away from the center toward higher and 
lower energies. These are identified in the semiclassical 
picture as the quantization of the electron motion along 
the equipotential line around the bottom or the top of 
the periodic potential. The widths of those levels be- 
come narrower for larger </> because the coupling between 
different potential valleys becomes exponentially small as 
the magnetic length becomes smaller. 

The Hall conductivity a^y is calculated by the Kubo 
formula as 



\vx\P)W\vy\a)-{a\vy\P){P\vx\a) 



(7) 

where Ep is the Fermi energy and the energy of the 
eigenstate \a) in the lowest Landau level. This is rewrit- 
ten in an ideal system as |2j 
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du„ 



dky 



, (8) 



where the summation is taken over all the occupied 
states. It is shown that the contribution from all the 
states in one subband is always quantized in units of 
—e^ /h 0. So we have an integer Hall conductivity when 
the Fermi energy Ep is in every gap, since it simply be- 
comes the summation of the integers for all the subbands 
below Ep. In such a gapped situation, we have a use- 
ful expression for the Hall conductivity called the Streda 
formula, which is directly derived from the Kubo 
formula Q. This is written as 



dn 

'as' 



(9) 



where n is the number of states below the Fermi energy 
per unit area and B is the magnetic field. In Fig. ^ 
we put the value of the Hall conductivity for each of the 
gaps obtained by this formula. The contribution by one 
subband, l^axy, is obtained as the difference in (Jxy be- 
tween the gaps above and below the subband. In Fig. 
(b) , the series of branches mentioned above always car- 
ries A(Ta;y = 0, since those subbands, coming from the 
localized orbitals around the potential minima or max- 
ima, contain a constant number of states n = 1 / o? which 
is independent of the magnetic field. 



III. LOCALIZATION AND THE QUANTUM 
HALL EFFECT 



We now move on to the disordered system to investi- 
gate the localization and the quantum Hall effect. Here 
we numerically diagonalize the Hamiltonian of finite 
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In the strong disorder regime, the energy scale of the 
disorder is characterized by 



FIG. 1: (a) Ideal energy spectrum of the lowest Landau 
level in a weakly modulated two-dimensional system, plotted 
against the magnetic flux (j). Integers inside gaps indicate the 
quantized Hall conductivity in units of ( — e^/Zi). (b) Zoom 
out of the figure (a). 



systems oi L x L with L ~ Ma (M is an integer), and 
calculate the localization length Lioc and the Hall con- 
ductivity axy as functions of Fermi energy. The localiza- 
tion length is obtained from the system size dependence 
of the Thouless number g{L), by assuming that g{L) be- 
haves as exp(— L/Lioc), where g{L) is defined as the ra- 
tio of the difference in the eigenenergies between periodic 
and antiperiodic boundary conditions to the mean level 
spacing AE. The Hall conductivity is calculated by the 
Kubo formula O with the exact eigenstates in the dis- 
ordered system, where the mixing between different Lan- 
dau levels is taken into account within the lowest order 
in (Vp -\- Vd)/{fiujc) with ujc — eB/m being the cyclotron 
frequency. Every quantity is averaged over a number of 
different samples. 



2ttP ' 



which represents the increasing width of the Landau level 
in absence of the periodic potential. |Q| This parameter 
is relevant as long as the disorder is strong enough to 
destroy the band structure inside the Landau level. In 
the weak disorder limit, on the other hand, the relevant 
parameter is the increasing width of the Bloch subband, 
which is written as 



SE w 2TrniVQp, 



(10) 



with p being the density of states per unit area. For 
a subband in the flux (j) = p/Qi P is approximated as 
\/{Wqa^) with the subband width W , so that we obtain 



5E 



2'iniiVQ 
Wqa? 



(11) 



The relation between T and 7 is given by 
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As a typical result, we show the calculation for the dis- 
ordered system in the flux = 3/2 in Fig. |21 The top 
panel indicates the system size dependence of Uxy and 
the density of states, and the bottom the inverse local- 
ization length l/Lioc. The parameter of the disorder is 
set to T /V — 0.25. In the clean limit, a Landau level 
splits into three separated subbands since p = 3, and 
the Hall conductivity carried by each of the subbands 
is AcTj-j, ~ 1,-1,1. Here and in the following we show 
Uxy and AcTxy in units of —e^/h. When the system is 
disordered we flnd that the gaps between the subbands 
are smeared out while the Hall conductivity converges to 
quantized values around the density of states (DOS) dips, 
indicating the appearance of the Hall plateau. If we look 
at the size dependence of the Hall conductivity in the 
whole energy region, we see that Uxy always approaches 
1 in the area Uxy > 1/2 in increasing the size, and in 
Uxy < 1/2, so we expect that in an infinite system the 
continuous function axy{E) changes into Hall plateaus 
connected by steps, as shown by the stcplikc line. There- 
fore we speculate that the points of cjxy — 1/2 are identi- 
fied as the extended levels in an infinite system, and they 
indeed agree with the energies where localization length 
diverges as shown in Fig. |21 (b). 

We can see similar tendencies in the size dependence 
of axy in other fluxes as well, where the flxed points are 
always found at axy ~ n+l/2 (n is an integer) in all the 
cases investigated. While each of the three subbands has 
one extended level at a certain energy in (j) — 3/2, the 
localization generally depends on the Hall conductivity 
assigned to subbands. We have shown that a subband 
carrying zero Hall conductivity is completely localized in 
finite disorder strength. Jljj 
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FIG. 2: (a) Hall conductivity a^y (solid line) and the density 
of states p (dashed) in the flux (p = 3/2 and the disorder 
r/V = 0.25 with various sizes L/a — 8, 18, and 36. The gray 
steplike line is an estimate of a^y as 1/ — > oo. (b) Inverse 
localization length estimated from the Thouless number. The 
vertical lines penetrating the panels represent the energies of 
the extended levels. 



IV. EVOLUTION OF THE EXTENDED LEVELS 

When we consider a situation where the disorder in- 
creases in the Hofstadter butterfly, we expect that the 
extended levels in the subbands move and merge on the 
energy axis in some way, and in large enough disorder, 
only one remains at the center of the Landau level. It 
would be nontrivial and intriguing to ask how they evolve 
as a function of disorder, and how different they are in 
various fluxes. Here we study the evolutions for some 
particular fluxes, on the assumption that the extended 
levels always exist at the energies of axy — n + 1/2 with 
integers n, and the region where 71 — 1/2 < axy < n+l/2 
in finite systems becomes a^y = n in an infinite system. 

We first show the results for the flux ^ = 3/2 in Fig. |31 
The top and middle panel indicate the density of states 
and the Hall conductivity, respectively, in a finite system 
L = 12a for several disorder parameters F. The bottom 
panel shows the traces of the extended levels obtained by 



just taking the energies of axy 



-1/2. We observe that 



three extended levels get closer as the disorder increases, 
and contract into one at a certain F. The combination 
of three branches at a time is due to the electron-hole 
symmetry with respect to E — 0, owing to the equal 
amounts of the positive and negative scatterers. If we 
break the symmetry by introducing imbalance between 
them, the evolution changes in such a way that two of 
them annihilate and one is left intact. |l5j | 

The situation becomes a little complicated in a slightly 
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FIG. 3: (a) Density of states and (b) the Hall conductivity as 
a function of the Fermi energy in = 3/2 and L = 12a with 
several disorder parameters T/V. (c) Trace of the extended 
levels as a function of disorder, obtained by taking the energies 
oiffxy = n-|-l/2 (in units of —e^/h). The horizontal bars show 
the numerical errors. Integers in the enclosed areas represent 
the quantized Hall conductivities. The thick bars at F = 
and the integers below them indicate the ideal subbands and 
Aa^y, respectively. 



larger flux — 5/3, as shown in Fig. ^ In the 
clean limit, we have five subbands carrying Aaxy = 
— 1, -1-2, — 1, -1-2, — 1 from the bottom to the top. Since 
the second and the fourth bands pass through n -I- 1/2 
twice, each of them comes to have two extended levels 
in sufficiently small disorder on the present assumptions. 
On going to stronger disorder, one of the two extended 
levels experiences a pair annihilation with the lowest or 
highest subband, and only three extended levels are left 
as in = 3/2. When we start with a more complicated 
fiux, we speculate that the Hall plateau structure is sim- 
plified one by one, going through analogs of the simpler 
fiuxes nearby. 

To check that there are actually two separated ex- 
tended levels in a subband carrying Aaxy = 2, we can 
make a scaling analysis in the Hall conductivity on vary- 
ing the system size. In the flux (j) ~ 5/3, unfortunately, 
the large statistical error prevents us from resolving the 
splitting, but we can in a similar situation in another flux 
(j) = 8/15, where the lowest subband carrying Aaxy — 2. 
Figure |31 shows the Hall conductivity and the localization 
length for the lowest subband with disorder F = 0.01731^. 
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FIG. 4: Plots similar to Fig. 01 calculated for = 5/3. 



The difference in axy between L/a ~ 45 and 30, shown 
in the middle panel, indicates that axy actually has turn- 
ing points around axy = 1/2 and 3/2. This is confirmed 
by the calculation of the localization length in the lowest 
panel, showing that Lioc diverges at those two energies, 
so it is quite likely that there are two extended levels in 
this subband in an infinite system. 

A major difference from = 5/3 is that the ideal sub- 
band width is much narrower than the energy gap in the 
present case, so we can set the disorder in such a way 
that the states in the subband are completely mixed up 
but not with other subbands. We speculate that this sit- 
uation makes the localization length smaller and enables 
us to resolve the separated extended levels in a finite size 
calculation. We have another example where one sub- 
band has more than one extended levels in an anisotropic 
modulated quantum Hall systems. There axy{E) in a 
subband with /S.axy = 1 has a non-monotonic behavior 
crossing axy = 1/2 three times, so that three extended 
levels arise. 

Last we investigate a flux = 8/7, slightly away from 
(j) — 1. Around this region the Landau level splits into 
a number of tiny subbands as seen in the original Hofs- 
tadter butterfly, just like the Landau levels in the usual 
2DEG around zero field. This is actually understood as 
the Landau quantization in the magnetic Bloch band at 
0=1 caused by the residual flux (j> — 1. As seen 
in Fig. El each of these "Landau levels" equally carries 
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FIG. 5: (a) Hall conductivity axy and the density of states 
p calculated for the lowest subband in = 8/15, carrying 
Aa^y = -2e'^/h, with the disorder T/V = 0.0173 and various 
sizes L/a = 30, 45, and 60. The gray steplike line is an esti- 
mate at L — > 00. The vertical line at i5 ~ —0.0865 represents 
the energy region of the ideal subband. (b) Relative value of 
axy in L/a = 60 estimated from L/a = 45. (c) Inverse lo- 
calization length (in units of a) estimated from the Thouless 
number. Two vertical lines penetrating the panels represent 
the energies of the extended levels in an infinite system. 



the Hall conductivity of Aaxy — 1 while the center one 
has a large negative value. It is highly nontrivial how the 
extended states evolves in increasing disorder in such a 
situation. 

Figure shows the disorder dependence of the density 
of states, the Hall conductivity axy, and the traces of 
<^xy = 71 + 1/2 in a fixed system size L = 14a in = 
8/7. Remarkably the result suggests that a pair creation 
of extended states can occur in increasing disorder, as 
indicated by the arrows in (c) . In (b) , we can see that this 
happens when a dip in axy (E) touches the line of axy = 
n + 1/2, where a pair of extended levels with opposite 
Hall conductivity are created as schematically shown in 
Fig. [71 In even stronger disorder, these newly created 
pairs annihilate with other extended levels in the lower 
and higher subbands. 

To confirm the existence of a pair creation, we investi- 
gate the scaling behavior of the Hall conductivity axy [E) 
in a similar situation in = 5/4, for the smaller disorder 
F = 0.4F and the larger 0.525F. The results in Fig. |H1 
show that the direction of the size dependence around 
the dip {E/V w —0.2) seems to change when passing 
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FIG. 6: Plots similar to Fig. |3|calculated for (j) — 8/7. Arrows 
indicate the positions of the pair creations of the extended 
levels. 
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FIG. 7: Schematic figure of pair creation of the extended 
levels. The left and right panels show o-xy{E) in smaller and 
larger disorder, respectively, where a solid curve is for a finite 
system and a gray line for an infinite system. The pair of 
extended levels are newly created when a dip touches the line 
of axy = 0.5(— e^//i). 



available. This is consistent with the proposed scenario 
as in Fig. [71 while we need further extensive work to 
confirm that the tendency remains towards the infinite 
system. 

For the usual 2DEG without a periodic potential, there 
has been a long debate about the fate of the extended 
states in the Landau levels in the strong disorder limit 
in a weak magnetic field. It has been argued that all the 
extended states float up in the energy axis with increas- 
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FIG. 8: Hall conductivity a^y in the fiux <^ = 5/4 and the dis- 
order r/V (a) 0.4 and (b) 0.525 (in different vertical scales), 
with the system sizes L/a — 16 and 36. 



ing disorder and the system finally becomes an insula- 
tor. [itI Us . 19.] Numerical calculations in tight-binding 
lattice models, on the other hand, showed that the the 
states with negative Hall conductivity move down from 
the tight-binding band center, to annihilate with the ex- 
tended states in the lower Landau levels. [20ll27. 22] Our 
theory may suggest a somewhat different story where pair 
creations of the extended states occur, since we actually 
observe similar dips in axy {E) around higher Landau lev- 
els in a disordered nonperiodic 2DEG. 

We also remark that the evolution of the extended lev- 
els generally depends on the correlation lengths of the 
disorder potential. It was shown in the tight-binding 
model that the pair annihilation of the extended states 
with positive and negative Hall conductivity do not occur 
in the long range disorder. |2^ In our problem, the 
scenario of the pair creation might change in a correlat- 
ing disorder, since the appearance of the d ips in axy{E) 
is due to the inter-Landau-level mixing [2^ , and thus 
is suppressed by the long-range scatterers. This should 
be addressed in future work. 



V. HALL PLATEAU DIAGRAM 

We extend the analysis to general fluxes to look into 
the quantum Hall effect over the whole region of the Hofs- 
tadter butterfly. Figure|Hlshows (a) the Hall conductivity, 
(b) the density of states, and (c) the localization length 
in disordered systems with T /V = 0.25, plotted against 
the Fermi energy and the magnetic field. The solid lines 
in Fig. 13 (a) represent the extended levels in an infinite 
system, and the enclosed areas are the Hall plateaus with 
the quantized value of axy^ where we again assume that 
the extended levels in an infinite system coincide with 
the energies of Uxy = n -I- 1/2 in a finite system. In the 
density of states [Fig. Elb)], we recognize valleys in the 
contour plot as the remnants of the minigaps, and some 
of them can be associated with the Hall plateaus in (a). 
For instance, a pair of valleys between <f> = \ and 2 inside 
the Landau level corresponds to two major plateaus with 
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cr^j^ = in <£; < 0.5V and a^y = 1 in -0.5V < E < 0. 

In Fig. I^c), we see that the energy region with large 
locahzation length becomes broad particularly around 
integer fluxes. There the cluster of subbands widely 
spreads along the energy axis, and thus the states are 
less easily mixed. We also find that the upper and lower 
branches of the spectrum in </> > 2, carrying zero net Hall 
conductivity, are all localized, as seen in the particular 
case (j) = 3. The localization length is smaller for 
larger (p in those branches, mainly because the mixing of 
the states becomes stronger in narrower subbands. 

In Fig. ^1 we show the Hall plateau diagrams and 
the density of states for the different disorder parameters 
T/V — 0.15 and 0.5, which exhibit the dependence on F 
together with Fig. Ufa) and^b). We can see that the 
small Hall plateaus coming from the fine gap structure 
gradually disappear as the disorder becomes larger, and 
the only extended level is left at the center of the Lan- 
dau level in the strong disorder limit. We here notice in 
T/V — 1.5 that the two largest plateaus between (p — 1 
and 2 mentioned above are more easily destroyed around 
(f) ~ 1 than around </> = 2, or the right end is detached 
from (p = 1 while the left sticks to = 2. This is because 
the tiny gaps around (p = 1 are easily swallowed up by 
the large density of states around the center. 



VI. SELF-SIMILARITY 

The interesting observation in the Hofstadter butterfly 
is that the identical spectrum with the identical distribu- 
tion of the Hall conductivity repeatedly appears in differ- 
ent fluxes. In Fig. ^b), we can see that the spectrum 
and the Hall conductivity at (p correspond to the middle 
part of 0-|-2 with the top and bottom branches excluded, 
so that the entire structure repeats over (p -\- 2m with in- 
teger m. The proof of the correspondence is presented in 
Appendix. We note that, while the Hofstadter butterfly 
has similar gap structures everywhere in a fractal fash- 
ion, the Hall conductivities in the corresponding gaps 
do not always coincide. An example of the incomplete 
correspondence is seen between <p and <p' = 1/(1— l/(p) 
(such as (p — i and 3/2). In the tight-binding model, it 
was shown that the distributions of the Hall conductivity 
within a cluster resemble each other up to a scale factor 
among some series of fluxes. 

In the following we show that, in the presence of dis- 
order, the corresponding clusters with identical Hall con- 
ductivity have qualitative agreement also in the localiza- 
tion length. Here we particularly take a pair of fluxes 
(p = 3/2 and 7/2, in both of which the central three 
subbands have the Hall conductivity (1,-1,1). In Fig. 
im we compare the disorder effects on the density of 
states and the localization length. We set the disorder as 
F = 0.25 and 0.20 for (p = 3/2 and 7/2, respectively, so 
that the renormalized DOS broadening 7/Wtot is equiv- 
alent, where we defined 7 by putting in Eq. the 
full width of three subbands Wtot • The result shows that 



r/V'=0,25 
(a) Hall conductivity 




(c) Localization iengtti 




FIG. 9: (a) Hall plateau diagram in the lowest Landau level 
at the disorder T/V — 0.25, plotted against the magnetic 
flux and the Fermi energy. Solid lines show the energy of 
(^xy = n -I- 1/2 (n is an integer), which are identified as the 
extended levels in an infinite system, and the integers indicate 
the Hall conductivity in units of —e^/h. The ideal spectrum 
is shown as the gray scale, (b) Corresponding plots for the 
density of states in units of l/{Va'^). (c) Localization length 
in units of a. The dashed lines indicate the extended levels 
shown in (a). 



the densities of states are broadened equivalently as ex- 
pected, and that the localization lengths iioc then agree 
qualitatively without any scale factors. 

Figure [TTT c'l shows the evolution of the extended lev- 
els as a function of the disorder, which are obtained by 
taking the points of axy = 1/2. Three extended levels 
in = 3/2 and 7/2 evolve in a parallel fashion with the 
disorder strength 7/ Wtot, where they come closer as 7 
becomes larger and combine into one at 7/Wtot ~ 0.4. 
The critical disorder at which the combination occurs is 
slightly smaller in (p — 7/2 than in 3/2, presumably be- 
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FIG. 10: Hall plateau diagram (left) and the density of states 
(right) in the disorder T/V — 0.15 (top) and 0.5 (bottom), 
plotted against the magnetic flux and the Fermi energy. The 
corresponding parameter region is shown as the dashed line 
in Fig. Ela). 



cause in = 7/2 the level repulsion from the outer sub- 
bands (out of the figure) pushes the states in the central 
three subbands toward the center of the spectrum and 
that enhances the contraction of the extended levels. We 
expect that the critical behavior of the three extended 
levels at the combing point is universal, but we could not 
estimate the critical exponent in this simulation due to 
statistical errors. The evolution of the plateau diagram 
should become basically similar among + 2m, so we 
know all from the information of the first unit. 

The similar subband structures with the identical Hall 
conductivity can be found in other hierarchies in the Hof- 
stadter butterfly. The smallest unit that we have found 
is schematically shown in Fig. 1121 where the structure 
indicated by A between = and 1 repeats over the 
entire region. We can show similarly that each unit has 
the identical distribution of the number of states to each 
subband so that we have the identical Hall conductivity 
for every gap. We expect that the localization length in 
each single unit A becomes similar when the disorder is 
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FIG. 11: (a) Density of states and (b) the inverse localization 
length calculated for disordered systems with </!> = 3/2 (solid 
lines) and 7/2 (dashed). The energy scale is normalized by the 
whole width of the Landau level in the ideal system, Wtot , and 
the density of states in units of l/(VFtota^). The parameter of 
the disorder is y/Wtot ~ 0.035. (c) Evolutions of the extended 
levels with changing disorder 7, in (/> = 3/2 (filled circles) and 
7/2 (triangles). 




FIG. 12: Schematic diagram showing decomposition of the 
Hofstadter butterfly into identical units, where the gap struc- 
ture and the Hall conductivity inside each gap coincide. The 
letter A indicates the direction. 



weak enough that the mixing among different units can 
be neglected. 
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VII. CONCLUSION 

We studied the quantum Hall effect in a Landau level 
in the presence of a two-dimensional periodic potential 
with short-range disorder potentials. It is found that, in 
all the cases we studied, the Hall conductivity becomes 
size independent at a^y = n + 1/2 (in units of —e^/h), 
and those points are identified as the extended levels in 
an infinite system. We propose a possible model for the 
evolution of the extended levels by tracing Uxy = n + 1/2, 
which predicts that a subband with axy = n has n (or 
more) bunches of extended levels, and possibly that pair 
creation of the extended levels can occur in increasing 
disorder, as well as pair annihilation. We also find that 
the clusters of subbands with an identical Hall conductiv- 
ity, which compose the Hofstadter butterfly in a fractal 
fashion, have a similar localization length in the presence 
of the disorder. 



ACKNOWLEDGMENTS 

This work has been supported in part by the 21st 
Century COE Program at Tokyo Tech "Nanometer- 
Scale Quantum Physics" and by a Grant-in- Aid for COE 
(12CE2004 "Control of Electrons by Quantum Dot Struc- 
tures and Its Application to Advanced Electronics" ) from 
the Ministry of Education, Science and Culture, Japan. 
One of the authors (M.K.) is supported by a Grant-in- 
Aid for Scientific Reseach from the Ministry of Educa- 



[1] D. R. Hofstadter, Phys. Rev. B 14, 2239 (1976). 

[2] D. J. Thouless, M. Kohmoto, M. P. Niglitingale, and M. 

den Nijs, Phys. Rev. Lett. 49, 405 (1982). 
[3] C. Albreclit, J. H. Smet, K. von Klitzing, D. Weiss, V. 

Umansky, and H. Schweizer, Pfiys. Rev. Lett. 86, 147 

(2001); Physica E 20, 143 (2003). 
[4] M. C. Geisler, J. H. Smet, V. Umansky, K. von Klitzing, 

B. Naundorf, R. Ketzmerick, and H. Schweizer, Phys. 

Rev. Lett. 92, 256801 (2004). 
[5] D. Pfannkuclie and R. R. Gerhardts, Phys. Rev. B 46, 

12606 (1992). 

[6] U. Wulf and A. H. MacDonald, Pliys. Rev. B 47, 6566 
(1993). 

[7] Y. Tan, J. Pliys. Condens. Matter 6, 7941 (1994); Phys. 

Rev. B 49, 1827 (1994). 
[8] Y. Hatsugai, K. Ishibashi, and Y. Merita, Phys. Rev. 

Lett. 83, 2246 (1999). 
[9] K. Yang and R.N. Bhatt, Phys. Rev. B 59, 8144 (1999). 
[10] B. Huckestein, Phys. Rev. Lett. 72, 1080 (1994). 
[11] B. Huckestein, Rev. Mod. Phys. 67, 357 (1995). 
[12] M. Koshino and T. Ando, J. Phys. See. Jpn. 73, (2004) 
3243. 

[13] P. Stfeda, J. Phys. C 15, L718 (1982). 



tion, Science and Culture, Japan. Numerical calcula- 
tions were performed in part using the facilities of the 
Supercomputer Center, Institute for Solid State Physics, 
University of Tokyo. 



APPENDIX: CONCIDENCE IN THE HALL 
CONDUCTIVITY 

The equivalence in the gap structure with the Hall con- 
ductivity between cj) and (j) + 2, discussed in Sec. IVII is 
explained as follows. The number of subbands in a Lan- 
dau level is given by the numerator of the magnetic flux 
(j), so that we have p bands in (f> — p/q and p+2q bands in 
(j> + 2. Each single subband in and (/) -f 2 consists of the 
equal number of states per unit area, l/{qa^), since two 
fluxes have a common denominator q so that they have 
equal foldings of the Brillouin zone. We can then see that 
each of the top and bottom branches in (j) + 2 contains q 
subbands, because each has a constant number of states 
as explained in Sec. II. Thus the number of sub- 
bands in the middle part in + 2 (the top and bottom 
branches removed) becomes (p + 2q) — 2 x q = p, which 
is equal to the total subbands in (j). Now we see that 
the corresponding spectra between (p and (j) + 2 have the 
same number of subbands with the the equal number of 
states, so we come to the conclusion that the Hall con- 
ductivity becomes identical between the corresponding 
gaps, by using the Streda formula ©. 



[14] T. Ando and Y. Uemura, J. Phys. Soc. Jpn. 36, 959 
(1974). 

[15] M. Koshino and T. Ando, Physica E, 29, 588 (2005). 
[16] M. Chang and Q. Niu, Phys. Rev. Lett. 75, 1348 (1995). 
[17] H. Levine, S. B Libby, and A. M. M. Pruisken, Phys. 

Rev. Lett. 51, 1915 (1983). 
[18] D. E. Khmelnitskii, Phys. Lett. 106A, 182 (1984). 
[19] R. B. Laughlin, Phys. Rev. Lett. 52, 2304 (1984). 
[20] T. Ando, Phys. Rev. B 40, 5325 (1989). 
[21] D. Z. Liu, X. C. Xie, and Q. Niu, Phys. Rev. Lett. 76, 

975 (1996). 

[22] D. N. Sheng and Z. Y. Weng, Phys. Rev. Lett. 78, 318 
(1997). 

[23] D. N. Sheng, Z. Y. Weng, and X. G. Wen, Phys. Rev. B 

64, 165317 (2001). 
[24] T. Koschny, H. Potempa, and L. Schweitzer, Phys. Rev. 

Lett. 86, 3863 (2001). 
[25] T. Ando, Y. Matsumoto, and Y. Uemura, J. Phys. Soc. 

Jpn. 39, 279 (1975). 
[26] T. Ando, J. Phys. Soc. Jpn. 55, 3199 (1986). 
[27] H. Ritter, Z. Phys. B 56, 185 (1984). 



